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We calculate the two-loop QCD corrections to the produc- 
tion of the neutral supersymmetric Higgs bosons via the gluon 
fusion mechanism at hadron colliders, including the contribu- 
tions of squark loops. To a good approximation, these ad- 
ditional contributions lead to the same QCD corrections as 
in the case where only top and bottom quark loops are taken 
into account. The QCD corrections are large and increase the 
Higgs production cross sections significantly. 



The search for Higgs particles is an important compo- 
nent of the experimental program at future high energy 
hadron colliders. As such it is vital to have reliable pre- 
dictions for the production rates both in the Standard 
Model (SM) and in the minimal supersymmetric exten- 
sion of the Standard Model (MSSM). The two-loop QCD 
corrections to the main production process, the gluon fu- 
sion mechanism Jjj, have been calculated in the SM in 
Refs. [|]J|] and later generalized to the quark contribu- 
tions in the MSSM ||,|. The corrections are large and 
positive, increasing the production rates significantly. 

The MSSM requires the introduction of two Higgs dou- 
blets leading, after spontaneous symmetry breaking, to 
two neutral CP-even (h and H), a neutral CP-odd (A) 
and two charged {H ± ) Higgs particles ||. While in the 
SM the dominant contribution to Higgs boson production 
in the gluon fusion mechanism originates from top and, to 
a lesser extent, bottom quark loops, in the MSSM there 
are additional contributions to the production of the CP- 
even Higgs bosons from scalar squark loops. These con- 
tributions can be neglected for very heavy squarks. How- 
ever, many supergravity-inspired models predict squark 
(in particular stop and sbottom squark) masses signifi- 
cantly below 1 TeV Q . In this case, squark loop contri- 
butions to the Higgs-gluon couplings can be of the same 
order, or even larger, as the standard quark contribu- 
tions, as was recently stressed in Refs. 0. 

In this letter, we present the O(a^) QCD corrections to 
the cross sections cr(pp ^>TL + X) of the fusion processes 
for the neutral CP-even Higgs particles H — h, H 
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Because of CP invariance, squark loops do not contribute 
to the production of the CP-odd Higgs boson in lowest 
order. The QCD corrections from squark loops are evalu- 



ated in the heavy squark limit, where the calculation can 
be simplified by extending the lowest-order low-energy 
theorems |||^] to two loops. This limit should be a very 
good approximation || for the production of Higgs par- 
ticles with masses smaller than twice the squark masses. 
Given the experimental bounds on the squark masses || , 
this is fully justified in the case of the lightest Higgs bo- 
son h, which is constrained to be lighter than ~ 130 GeV 
in the MSSM; for the heavier CP-even Higgs boson H, 
this approximation is valid for masses smaller than a few 
hundred GeV. For simplicity, we will restrict ourselves to 
the case of degenerate squarks where mixing effects are 
absent (in the absence of gluino-exchange, the results can 
be trivially generalized to include mixing). Also in this 
case, scalar squarks will not contribute to the production 
of the CP-odd Higgs boson A at next-to-leading order. 

To lowest order, the cross sections for CP-even Higgs 
boson production at proton colliders are given by 
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with dC" /dr-H the gluon luminosity at t-h — Mh/s and 
s is the total cm. energy. The parton cross sections are 
built up from heavy quark and squark amplitudes, 
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where the sums run over t, b quarks and the left- and 
right-handed squarks Ql,Qr, which in the absence of 
mixing are identical to the mass eigenstates. The form 
factors, with the scaling variables 7q/q = Am^^/M^, 
can be expressed as 



MTQ)=T Q [l+(l-T Q )f(T Q )} 
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The normalized scalar quark and squark couplings to the 
CP-even Higgs bosons, =, can be found in Refs. UK- 
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In the case where all squarks are taken to be degenerate, 
only the contributions proportional to the Yukawa-type 
couplings of the stop and sbottom squarks have to be 
added to the top and bottom quark loop contributions. 
The couplings as well as the CP-even Higgs masses are 
determined at tree-level by two parameters, which are 
generally chosen to be the ratio of the vacuum expecta- 
tion values of the two Higgs fields, tg/3, and the pseu- 
doscalar Higgs mass, Ma- All MSSM Higgs masses and 
couplings are calculated using the two-loop renormaliza- 
tion group improved effective potential |l0[ . 

The QCD corrections to the gluon fusion process, 
eq.(Q), consist of virtual two-loop corrections and one- 
loop real corrections due to gluon radiation, as well as 
contributions from quark-gluon initial states and quark- 
antiquark annihilation. The renormalization program 
has been carried out in the MS scheme for the strong cou- 
pling constant and the parton densities, while the quark 
and squark masses arc defined at the poles of their re- 
spective propagators. The result for the cross sections 
can be cast into the form: 



a(pp ->H + X) = a, 
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The coefficient C~n denotes the virtual two-loop correc- 
tions regularized by the infrared singularities of the real 
gluon emission. The terms Act^ = g.q) denote the 
finite parts of the real corrections due to gluon radiation 
and the gq and qq initial states. The expressions for the 
t, b quark contribution can be found in Refs.||-Q]. 

The calculation of the QCD corrections has been per- 
formed by extending the low-energy theorems |||| to 
scalar squarks at the two-loop level. For a light CP-even 
Higgs boson, these theorems relate the matrix elements of 
the quark and squark contributions to the Higgs-gluon 
vertex to the gluon two-point function. Denoting the 
matrix element of the squark contribution to the gluon 
two-point function by Mg(gg) and the corresponding 
matrix elements with an additional light CP-even Higgs 
boson by MQ{ggT(), one has at lowest order |flj 
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To extend this relation to higher orders, all quantities 
have to be replaced by their bare values; after differenti- 
ation, the renormalization then has to be performed. In 
the following we consider only the pure gluon exchange 
contributions, which are expected to be the dominant 
ones; for heavy enough gluinos, the two-loop corrections 
due to gluino exchange should be small since they are 
suppressed by inverse powers of the gluino mass. In this 
case the differentiation with respect to the bare squark 
mass m% can be rewritten in terms of the renormalized 



mass m,Q. A finite contribution to the QCD corrections 
arises from the anomalous mass dimension ]3| 
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The remaining differentiation with respect to the renor- 
malized squark mass of the gluon two-point function 
leads to the squark contribution /3q to the QCD f3 func- 
tion. The final result for the squark contributions to the 
H coupling to gluons can be expressed in terms of the 
effective Lagrangian 
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The QCD corrections are then fully determined by the 
anomalous mass dimension of the squarks PJl^l 
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and the squark contribution to the QCD (3 function (U 
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resulting in a final rescaling of the lowest-order La- 
grangian by a factor 1 + 25a s /67r at next-to-leading 
order, compared to a rescaling factor 1 + lla s /47r for 
the quark contribution. Starting from the Lagrangian 
eq.(10), the effective QCD corrections due to real gluon 
emission and the gq/qq initial states have to be added. 
These corrections are identical to the corresponding cor- 
rections to quark loops || in the heavy quark limit. 

The QCD corrected squark loop amplitudes have to be 
added coherently to the corrected t, b loop amplitudes, 
whose full mass dependence is known. To obtain a more 
reliable prediction for the total cross sections, the result- 
ing amplitudes for the squark contributions have been 
normalized to the lowest-order amplitude in the limit of 
large squark masses. These ratios are then multiplied 
by the lowest-order amplitude including the full squark 
mass dependence. The heavy squark limit is then ex- 
pected to be a very good approximation for Higgs masses 
below the QQ* threshold, as in the corresponding case of 
top quark contributions 

The final results for the partonic cross sections defined 
in eq.(R) can be expressed as 
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with f = t-j-i/t and JV^ being the number of light fla- 
vors contributing to the evolution of a s and the parton 
densities. The renormalization scale /i enters the lowest- 
order expression as the scale of the strong coupling 
a s = a s (/i 2 ). Pag,P gq denote the Altarelli-Parisi split- 



ting functions |l4j ; M is the factorization scale at which 
the parton luminosities are evaluated. F + is the usual + 
distribution, F(f) + = F(f) - 5(1 - f) dxF(x). 

The contributions to the coefficients cY and d 1 ^ ap- 
pearing in eq.(13) from squarks, in the heavy squark limit 
without t, b loops, are given by 
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FIG. 1. K factors of the cross sections a(pp — > TL + X) for 
tg/3 = 1.5 and 30. The solid lines include t, b as well as squark 
contributions, the dashed lines include only the t, b contribu- 
tions. The common squark mass is chosen to be = 200 
GeV. We take mt — 5 GeV, mt = 176 GeV and use the 
next-to leading order a s , fixed by the world average value 
a s (M|) = 0.118 @. The cross sections are convoluted with 
next-to-leading order GRV parton densities [Q. The renor- 
malization scale [i and the factorization scale M are identified 
with the Higgs masses. 



In Fig.l, we present the K factors for the QCD cor- 
rections to the production of the CP-even MSSM Higgs 
bosons as functions of the Ti masses for the LHC at a 
cm. energy ^fs = 14 TeV with (solid lines) and with- 
out (dashed lines) the squark contributions. The K fac- 
tors are defined as the ratios of the QCD corrected and 
lowest-order cross sections, using next-to-leading order 
a s and parton densities in both terms. A common value 
mQ = 200 GeV has been used for the left- and right- 
handed stop and sbottom squark masses. This value is 
identified with the SUSY scale of the MSSM couplings 
and Higgs masses, leading to a rather low upper limit on 
the lightest Higgs mass for a given value of tg/3. 

The QCD corrections enhance the cross sections by 
a factor between 1.6 and 2.8, if the lowest-order cross 
sections are evolved with next-to-leading order a s and 
parton densities. If the lowest-order cross sections are 
convoluted with lowest order a s and parton densities, the 
K factors are reduced to a level between 1 and 2. It can 
be inferred from Fig.l that the inclusion of squark loops 
in the production of both CP-even Higgs particles h and 
H does not substantially modify the K factors compared 
to the case where squark loops are absent. The O(10%) 
discrepancy between the two factors for small tg/3 (where 
both the top and the stop contributions are dominant 
and can be approximated by their heavy mass limits) is 
mainly due to the difference between the contribution of 
quarks and squarks to the effective Lagrangian (c^ = 
for Q and = ^ for Q loops). 

We have verified that the K factors do not depend sig- 
nificantly on the squark mass which enters the MSSM 
couplings and lowest-order cross sections in our analy- 
sis. In fact, in the extreme situation where one of the 
t, b squark cigenstates is relatively light while the other 
squarks are heavy and decouple (as is the case for large 
squark mixing), the K factors are almost the same as 
in Fig.l. Therefore, while it substantially changes the 
Higgs-squark couplings and hence the production cross 
sections, mixing in the stop or sbottom sectors should 
have a rather modest impact on the K factors. 

Thus, to a good approximation, the effect of the squark 
loops in the gluon fusion mechanism is quantitatively 
determined by the lowest-order cross section (including 
squark loop contributions), multiplied by the known K 
factors when only the t, b quark contributions are 
taken into account. 

In Fig. 2, we illustrate the effect of including the squark 
loops and the QCD corrections to the production rate of 
the lightest CP-even Higgs particle. The ratio of the 
QCD corrected cross sections with and without squark 
loops is shown as a function of the common squark mass 
for three values of tg/3 = 1.5,3 and 30, with the pseu- 
doscalar Higgs mass fixed to Ma = 100 GeV. The squark 
contributions increase the production cross sections sig- 
nificantly for squark masses below about 500 GeV, espe- 
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cially for small and moderate values of tg/3; for higher 
masses tUq, the squarks decouple from the amplitude. 
This large effect can be understood by recalling that for 
these values of tg/3, the top and stop contributions dom- 
inate and the inclusion of the stop loops leads, in the 
heavy top and squark limit, to the enhancement of the 
lowest-order production cross section by an amount 



J t+t 
ft 




(15) 



where f t _i_j (ft) denotes the cross section including (with- 
out) stop loops. Thus for stop masses of the order 
of the top mass mt the cross sections are significantly 
enhanced by including stop loops. For large values of 
tg/3, because the Higgs couplings to (s)top (s)quarks are 
strongly suppressed (except for h in the decoupling limit), 
and the contribution of the bottom loop is enhanced by 
large logarithms compared to the sbottom loop, the pro- 
duction cross section is significantly affected only by light 
squark contributions, while they become negligible for 
sbottom masses of the order of 200 GeV. 
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FIG. 2. Ratio of the QCD corrected 
a(pp — > h + X) with and without squark loops for three values 
of tg/3 = 1.5,3,30, and for M A = 100 GeV. The secondary 
axes present the corresponding Higgs masses Mh. The quark 
masses, a s and the parton densities are as in Fig.l. 

Light squarks, with masses rrig < 300 GeV, can have 
a large impact on the production cross sections of the 
CP-even MSSM Higgs bosons at hadron colliders. We 
have presented the QCD corrections to the gluon fusion 
processes pp — > gg — > h/H including squark loops, in the 
heavy squark limit, which should well describe the full 
corrections including mass effects, at least in the range 
where the Higgs masses are smaller than twice the squark 
masses. The t and b quark mass dependence has been 
included exactly. To a good approximation the K factors 
are the same as the corresponding K factors when only 
top and bottom contributions are included. Since they 



increase the production cross sections significantly, the 
QCD corrections must be taken into account. 
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